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Abstract. In this article, we generalize the theory of discrete La- 
grangian mechanics and variational integrators in two principal directions. 
First, we show that Lagrangian submanifolds of symplectic groupoids give 
rise to discrete dynamical systems, and we study the properties of these 
systems, including their regularity and reversibility, from the perspective 
of symplectic and Poisson geometry. Next, we use this framework — along 
with a generalized notion of generating function due to Tulczyjew — to 
develop a theory of discrete constrained Lagrangian mechanics. This 
allows for systems with arbitrary constraints, including those which are 
"nonholonomic" (in an appropriate discrete, variational sense). In addi- 
tion to characterizing the dynamics of these constrained systems, we also 
develop a theory of reduction and Noether symmetries, and study the 
relationship between the dynamics and variational principles. Finally, we 
apply this theory to discretize several concrete examples of constrained 
systems in mechanics and optimal control. 



1. Introduction 

Understand [ing] in depth the symplectic geometry 
of the generating function method . . . should lead 
to a method for nonholonomic constraints. 

McLachlan and Scovel, A survey of open problems 
in symplectic integration jl6j 

In this article, we generalize the theory of discrete Lagrangian mechanics 
and variational integrators in two principal directions. First, we show that 
Lagrangian submanifolds of symplectic groupoids give rise to discrete dynam- 
ical systems, and we study the properties of these systems, including their 
regularity and reversibility, from the perspective of symplectic and Poisson 
geometry. Next, we use this framework — along with a generalized notion 
of generating function due to Tulczyjew — to develop a theory of discrete 
constrained Lagrangian mechanics. This allows for systems with arbitrary 
constraints, including those which are "nonholonomic" (in an appropriate 
discrete, variational sense). In addition to characterizing the dynamics of 
these constrained systems, we also develop a theory of reduction and Noether 
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symmetries, and study the relationship between the dynamics and varia- 
tional principles. Finally, we apply this theory to discretize several concrete 
examples of constrained systems in mechanics and optimal control. 

Before giving a more in-depth overview of the paper, and summarizing 
the key results, let us first give a brief background on discrete Lagrangian 
mechanics, in order to provide the context for the present work. 

1.1. Background. The subject of discrete Lagrangian mechanics concerns 
the study of certain discrete dynamical systems on manifolds. As the name 
suggests, these discrete systems exhibit many geometric features which are 
analogous to those in continuous Lagrangian mechanics: in particular, the 
dynamics of these systems satisfy variational principles, have symplectic or 
Poisson flow maps, conserve momentum maps associated to Noether-type 
symmetries, and admit a theory of reduction. While discrete Lagrangian 
systems are quite mathematically interesting, in their own right, they also 
have important applications to structure-preserving numerical simulation of 
dynamical systems in geometric mechanics and optimal control theory. 

In the simplest form of discrete Lagrangian mechanics, one begins with 
a function L: Q x Q ^ called the discrete Lagrangian, where Q is some 
configuration manifold. A trajectory qQ,...,qn G Q is a solution of the 
system if it extremizes the discrete action sum, 



where the endpoints go and qn are held fixed. This is essentially a discrete 
version of Hamilton's principle of stationary action, where the tangent bundle 
TQ has been replaced hy Q x Q, and the action integral has been replaced 
by an action sum. Solving this variational principle, one obtains the discrete 
Euler-Lagrange equations. 



for k = l,...,n— 1, where OqL and diL denote the partial derivatives 
of L with respect to the first and second arguments, respectively. These 
discrete Euler-Lagrange equations define an implicit function (i.e., a relation), 
mapping {qk-i,Qk) (,QkiQk+i)- If L is sufficiently nondegenerate, then this 
relation is the graph of a flow map QxQ^QxQ,at least locally. 

Alternatively, but equivalently, one can define a pair of discrete Legendre 
transformations F^L: Q x Q ^T*Q, given by 



F L{qk-i,qk) = -doL{qk-i,qk) , ¥+L{qk-i,qk) = diL{qk-i,qk) ■ 
This defines a symplectic relation on T*Q, given by 



n 




k=l 



= -S-S (go, • ■ ■ 
oqk 



qn) = doL {qk, qk+i) + diL {qk-i,qk) , 



Pk-i=¥ L{qk-i,qk), Pk =¥+L {qk-i,qk) 
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SO if the discrete Legendre transforms are (locally) invertible, then this gives 
a (local) symplectic automorphism 

F+L o {F-L)-' : T*Q ^ T*Q, (<7fc_i,Pfe-i) ^ {qk,Pk) , 

which can be interpreted as the discrete Hamiltonian flow of the system. 
More precisely, if lo denotes the canonical symplectic form on T*Q, then 
L:QxQ— s-Misa generating function for the Lagrangian submanifold 
dL {Q X Q) C {T*Q, —u) x {T*Q,ui), which can be identified with the graph 
of the discrete Hamiltonian flow in the nondegenerate case. Observe that 

F+L(gfc_i,gfc) =Pk= ¥'L{qk,qk+i) <^ doL{qk,qk+i) + diL {qk-i, qk) = 0, 

so this approach is equivalent to the previous one. 

Numerical methods which are constructed in this fashion are called vari- 
ational integrators, due to the key role played by the variational principle. 
This approach to discretizing Lagrangian systems was put forward in sem- 
inal papers by Suris [l9], Moser and Veselov [l7j, and others in the early 
1990s, and the general theory was developed over the subsequent decade (see 
Marsden and West [T^ for a comprehensive overview). 

Weinstein [22j observed that these systems could be understood as a special 
case of a more general theory, describing discrete Lagrangian mechanics on 
arbitrary Lie groupoids. Given a Lie groupoid G ^ Q and a discrete 
Lagrangian L : G — )■ M, a composable sequence of elements gi, gn & G is 
a solution trajectory if it extremizes the discrete action sum 

n 

5 {gi, ...,£/„) = ^L(c/fe) , 

k=l 

where the product gi • • • gn = 9 is held fixed. This variational principle yields 
a discrete Lagrangian flow map gk ^ gk+i- Equivalently, L can be interpreted 
as a generating function for the Lagrangian submanifold dL(G) C T*G. Since 
T*G is a symplectic groupoid over the dual Lie algebroid A*G, which has a 
canonical Poisson structure, this Lagrangian submanifold defines a Poisson 
relation on A*G. In the special case of the so-called pair groupoid QxQ ^ Q, 
this recovers the previous theory. We will recall the details of this approach 
in Section 2 , along with some of the more recent developments that appeared 
in Marrero et al. [I3j, including an explicit characterization of the discrete 
Lagrangian and Hamiltonian dynamics. 



1.2. Organization of the paper. In Section 2 we begin by recalling the 
basic definitions and properties of Lie groupoids and algebroids, highlighting 
the close relationship between symplectic groupoids and Poisson manifolds. 
After laying these foundations, we next describe how Lagrangian submanifolds 
of symplectic groupoids give rise to discrete dynamics, and how this can be 
seen as an abstract generalization of discrete Lagrangian mechanics. The 



main result of this section. Theorem 2.15, is a significant generalization 



of earlier results characterizing the regularity and reversibility of discrete 
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Lagrangian systems. We prove, in a generalized sense, that if either of the 
discrete Legendre transforms is a (local) diffeomorphism, then both are, and 
hence the discrete flow is a (local) Poisson automorphism. 



Next, in Section 3, we develop a theory of discrete constrained Lagrangian 



mechanics, for systems which are restricted to some constraint submanifold of 
a Lie groupoid. In this case, rather than having a discrete Lagrangian defined 
on the entire groupoid, we only have L: N ^ M for some N G G. When N is 
a proper submanifold, this cannot be a generating function in the usual sense; 
however, we show that a more general notion of generating function, due to 
Tulczyjew, can be used to generate a Lagrangian submanifold S/, C T*G. 
Using this approach, we characterize the regularity of discrete constrained 
Lagrangian systems, and obtain explicit formulas for the discrete Legendre 
transforms and discrete Euler-Lagrange equations. Finally, after showing 
that 'El C T*G has the structure of an affine bundle over N associated to 
the conormal bundle i'*N, we discuss how the dynamics on S/, can thus be 
interpreted as the evolution of configurations (on the base) and Lagrange 
multipliers (on the fibers). 

In Section 4} we introduce morphisms of discrete constrained Lagrangian 
systems, and use this to develop the reduction theory for these systems. 
The notion of Noether symmetries is also introduced, and we prove the 
appropriate discrete version of Noether's theorem, relating symmetries to 
conserved quantities. 

Subsequently, in Section 5 we show how the dynamics of discrete con- 
strained Lagrangian systems can be derived from a constrained variational 
principle on N, and demonstrate that these dynamics are equivalent to those 
obtained previously via the generating function approach. 

After this, in Section 6, we examine examples of discrete constrained 
Lagrangian systems for several applications. These applications include 
constrained mechanics and optimal control on Lie groups, an optimal control 
problem with nonholonomic constraints known as the plate-ball system, 
and time-dependent constrained mechanics with either fixed or adaptive 
time-stepping. 

Throughout the paper, we have occasion to draw on certain technical 
results involving bisections of Lie groupoids, which are also instrumental 
to the multiplicative structure of the cotangent groupoid and to the proof 
oflTheorem 5.11 We provide some supplementary details and discussion of 



these results in 



Appendix A 



2. SyMPLECTIC groupoids AND DISCRETE POISSON DYNAMICS 

2.1. Lie groupoids and algebroids. Let us first briefly recall some deflni- 
tions pertaining to Lie groupoids and Lie algebroids. This will provide an 
abbreviated reference and flx the notation used throughout the remainder of 
the paper. For a detailed treatment of this rich topic, see the comprehensive 
work of Mackenzie [11] . 
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Definition 2.1. A groupoid is a small category in which every morphism is 
invertible. That is, the groupoid denoted G ^ Q consists of a set of objects 
Q, a set of morphisms G, and the following structural maps: 

(i) a source map a: G ^ Q and target map {3 : G ^ Q; 

(ii) a multiplication map m : G2 — )■ G, {g, h) 1— )■ gh, where 

G2 = G pX^G = {{g,h)&GxG\P{g) = a{h)} 

is called the set of composable pairs, such that multiplication is 
associative whenever defined; 

(iii) an identity section e: Q ^ G oi a and /3, such that for all g E G, 

^{oi{g))g = g = ge {l3{g)) ; 

(iv) and an inversion map i: G G, g ^ g~^, such that for all g E G, 

gg-^ = e{a{g)), 9~^g = e{f3{g)) . 

Remark 2.2. Intuitively, groupoids can be viewed as either "categories with 
extra structure" or as "groups with missing structure," and both perspectives 
are useful. For instance, it can be helpful to see G G as an "arrow" from 
a{g) to P{g), where the multiplicative structure defines the composition of 

arrows • ""F^ • ""^i • . Alternatively, the groupoid G can be thought of 
as being a weaker version of a group, where multiplication is defined only 
partially (on G2 C G x G) rather than totally (on G x G). 

Next, for any element (7 G G of a groupoid, there are associated left and 
right translation maps, which act, respectively, on those elements with which 
g is composable on the left and on the right. 

Definition 2.3. Given a groupoid G ^ Q and an element g e G, de- 
fine the left translation £g: (l^id)) — ^ {(^{9)) right translation 
rg : {a{g)) ^ (Pig)) by g to be 

lg{h) = gh, rg{h) = hg. 

For discrete mechanics, we will focus on a particular class of groupoids. Lie 
groupoids, which — much like Lie groups, of which they are a generalization — 
have a differential structure in addition to (and compatible with) their 
algebraic structiure. 

Definition 2.4. A Lie groupoid is a groupoid G ^ Q where G and Q are 
differentiable manifolds, a and (3 are submersions, and the multiplication 
map m is differentiable. 

Remark 2.5. Observe that it is necessary for a and /3 to be submersions so 
that G2 is a differentiable manifold; otherwise, one could not make sense of 
m as a differentiable map. It also follows from the definition that m is a 
submersion, e is an immersion, and i is a diffcomorphism. 

In a natural way, one can now introduce the notion of left- and right- 
invariant vector fields on a Lie groupoid. 
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Definition 2.6. Given a Lie groupoid G ^ Q, a vector field ^ G X(G) is 
left-invariant if it is a-vertical (i.e., Ta{^) = 0) and Thig {^{h)) = ^ {gh) 
for all {g,h) G G2. Similarly, ^ is right-invariant if it is /9-vertical (i.e., 
rm = 0), and nrg (^(/i)) = ? (hg) for all (/i,^) G G2. 

Just as the "infinitesimal version" of a Lie group is a Lie algebra, the 
infinitesimal version of a Lie groupoid is a Lie algebroid. We will first define 
Lie algebroids as abstract structures, in their own right, before subsequently 
recalling the definition of the Lie algebroid associated to a particular Lie 
groupoid. 

Definition 2.7. A Lie algebroid is a real vector bundle A ^ Q, equipped 
with a Lie bracket •] on its space of sections r(^), and with a bundle 
map p: A ^ TQ called the anchor map. Furthermore, if we also denote by 

p: T{A) — > X{Q) the homomorphism of (7°° ((5)-modules of sections induced 
by the anchor, then we require that this satisfies the "Leibniz rule," 

lX,fYj=flX,Y}+p{X) [f]Y, 

for all X,Y e r{A) and / G (7°°(Q). 

Remark 2.8. Using the Jacobi identity, it is not hard to show that the 
induced map p: r(A) X{Q) is a Lie algebra homomorphism, where X{Q) 
is endowed with the usual Jacobi-Lie bracket [•, •] for vector fields. 

Definition 2.9. Given a Lie groupoid G ^ Q, the associated Lie algebroid 
AG Q is defined by its fibers AqG = kerr^(q)a, i.e., the space of a- vertical 
tangent vectors at the identity section, for each q £ Q. There is a bijection 
between sections X G F {AG) and left-invariant vector fields ^ G X{G), 
defined by 

(1) ^(5) = r,(^(,))^,(x(/?(5))). 

The Lie algebroid structure on AG is then given by the bracket [•, •] and 
anchor p satisfying the conditions 

ix,Yj = [XM, P{X)iq) = r,(,)^(X(g)), 
for all X,Y eF {AG) and qeQ. 

Remark 2.10. Alternatively, one can also establish a bijection between sections 
X G F {AG) and right-invariant vector fields ^ G X{Q), defined by 

(2) ^(5) = -r,(,(,))(r,oz)(X(a(5))), 
which yields the Lie bracket relation 

lX,Y\ = -[jl,f]. 

Thus, X I— )• ^ is a Lie algebra isomorphism, while X 1— )• ^ is a Lie algebra 
anti-isomorphism. 
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Note that for every v £ AqG, we have the following expression for the 
tangent to the inversion map: 

(3) T,^g)i{v) = -v + T,^g){eo f3){v). 

This can be seen by taking v = g{0) for some a- vertical path g{t), and observ- 
ing that T(^e{q),e{q))m{0,v) = {e{q) , g{t))\t=o = g{0) = v; likewise, for the 
/3-vertical path g~^{t), we have T(^e{q),e{q))iT^ O) = T^(^q^i{v). There- 

fore, T^e{q),e{q))fn {T^(^g)i{v) , v) =v + T^(^g)i{v)—hut since m {g-'^{t), g{t)) = 
€ {13(g)), we obtain v + Tg(q)i(t;) = T^(g) (e o /?) (v), which rearranges to give 
the desired equality. 

Finally, we discuss a subclass of Lie groupoids having even more additional 
structure. These are the symplectic groupoids, which are endowed with a 
symplectic manifold structure that is "compatible" with the Lie groupoid 
structure, in a sense which we will now define precisely. 

Definition 2.11. A symplectic groupoid is a Lie groupoid G ^ P, such that 

(i) (G, uj) is a symplectic manifold, 

(ii) the graph of m : G2 — )• G is a Lagrangian submanifold of G~ x G^ x G, 
where G~ = {G, —u) denotes the negative symplectic structure. 

If G ^ P is a symplectic groupoid and cu is the symplectic 2-form on G, 
then one may prove that (kerT^S)^ = kerT^/? for fj, £ G, where (ker T^a)'^ 
denotes the symplectic orthogonal complement of the subspace ker T^a. That 
is, the subspaces of a-vertical and /3-vertical tangent vectors are symplectic 
orthogonal to one another. Moreover, there exists a unique Poisson structure 
on P such that 5: G — )• P is anti-Poisson and ^5: G — )• P is Poisson. In 
addition, the identity section P — )■ G is a Lagrangian immersion. (See 
Coste et al. [4J, Mackenzie [11], Marie p^.) 

One example of a symplectic groupoid, in particular, lies at the heart of 
discrete Lagrangian mechanics: this is the cotangent groupoid. Given a Lie 
groupoid G ^ Q, let A*G ^ Q he the dual vector bundle of the associated 
Lie algebroid AG; then the cotangent groupoid is a symplectic groupoid 
T*G ^ A*G. Given an element n G T*G, the source and target are defined 
such that, for all sections X £T (AG), 

(4) {a{f,),X{aig))) = {fi,l{g)), {M,X{^{g))) = {fi,1({g)). 
(The definition of multiplication in r*G is slightly more intricate, and involves 



left and right translation by bisections; these are discussed in Appendix A 



2.2. Discrete dynamics of Lagrangian submanifolds. In this section, 
we discuss how a Lagrangian submanifold S C G of a symplectic groupoid 
G ^ P gives rise to discrete dynamics. In general, the dynamics are only 
defined implicitly, as a relation, rather than as an explicit flow map; we 
describe the conditions under which an explicit dynamical flow can be given. 
These dynamics can be interpreted either as discrete Lagrangian dynamics 
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on S or as discrete Hamiltonian dynamics on P. In the discrete Hamiltonian 
case, the exphcit flow is shown to be a (local) Poisson automorphism. 

Definition 2.12. Given a symplectic groupoid G ^ P, let S C G be 
a Lagrangian submanifold. Then a sequence ni, . . . ,fin £ G satisfies the 
discrete Lagrangian dynamics of S if /Ui, . . . , /u„ G S and 

^ (fik) = a (nk+i) , k = l,...,n-l, 

i.e., /Ui, . . . , /i„ E S forms a composable sequence in G. 

Example 2.13. Let G =^ Q be a Lie groupoid, equipped with a function 
L : G — 7- M called the discrete Lagrangian. Then L generates a Lagrangian 
submanifold dL(G) C T*G of the cotangent groupoid T*G ^ A*G. A 
sequence ^i, . . . £T*G satisfies the associated Lagrangian dynamics if 
and only if 

fik = dL (gk) for some gk ^ G, k = 1,. . .,n, 

and _ 

(/3 o dL) (gk) = (5 o dL) (gk+i) , k = l,...,n-l. 

Applying the definition of 5,/3 on T*G, this means that for any section 
X gT (AG), this is equivalent to 

{dL (gk) ,^ (gk)) = (dL (5^+1) , ^ A; = 1, . . . , n - 1, 

which can also be written as 

^[L] (gk) = 1[L] (gk+i) , A; = 1, . . . , n - 1. 

These are precisely the discrete Euler-Lagrange equations obtained by Mar- 
rero et al. [T^ for discrete Lagrangian mechanics on Lie groupoids. In the 
special case where G is the pair groupoid Q x Q ^ Q, this recovers the 
discrete Euler-Lagrange equations, 

doL {qk, qk+i) + diL {qk-i,qk) =0, k = 1, . . . ,n - 1, 

as in Marsden and West [TSl. 



These dynamics are implicitly defined, since they are given by a relati on 



Definition 2.12 



m 



on G rather than an explicitly defined map. Restating 
these terms, we see that m, . . . ,fin E G satisfies the discrete Lagrangian 
dynamics if and only if each pair of successive elements satisfies the relation 
i/J-k, IJ-k+i) € G2 n (S X S). (In fact, such a relation may be defined by any 
subset of a groupoid, although in general, the resulting dynamics will not 
preserve any notable geometric structure.) 

This raises the following question: Under what conditions is this relation, 
in fact, the graph of an explicit fiow map fik '-^ fJ-k+i (at least locally), and 
what properties does this fiow map have? Clearly, if the restricted source 
map 5|s: S — )• P is a (local) diffeomorphism, then this fiow is given by 
the composition (a|s)~^ o Furthermore, if the restricted target map 
S — )• P is also a (local) diff'eomorphism, then the flow is reversible. 
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and its (local) inverse is (/^Is) ""^ o ajs- In case both the restricted source 
and target maps are (local) diffeomorphisms, we say that S is a (local) 
Lagrangian bisection of the symplectic groupoid G, since it is both a (loc al) 



for 



a- and /3-section, as well as a Lagrangian submanifold. (See Appendix A 
more on bisections of Lie groupoids.) 

We now show that, in fact, if either of the maps a\-£ or /3|s is a local 
diffeomorphism, then both are. That is, if the Lagrangian submanifold T, C G 
is either a local a- or /3-section, then it is both, i.e., S is a local Lagrangian 
bisection. (Note that this result depends crucially upon the fact that G 
is a symplectic groupoid and that S C G is Lagrangian; it is not true for 
arbitrary submanifolds of Lie groupoids.) We begin by proving a lemma on 
complementary subspaces in a symplectic vector space; in the main theorem, 
this lemma is then applied to the a- and /3-vertical tangent subspaces, 
exploiting the fact that they are symplectic orthogonal complements. 

We note that this is a significant generalization of previous results on 
the regularity of discrete Lagrangian dynamics. Weinstein [22] first raised 
the question of how regularity results for the pair groupoid Q x Q might 
be generalized to arbitrary Lie groupoids G ^ Q, and this question was 
answered by Marrero et al. |13^ Theorem 4.13]. Here, we extend this answer 
from the Lagrangian submanifold dL(G) C T*G to arbitrary Lagrangian 
submanifolds of symplectic groupoids. 

Lemma 2.14. Let {Z,uj) be a Id-dimensional symplectic vector space with 
subspaces V,W C Z , and denote their symplectic orthogonal complements 
by V^, C Z, respectively. If dim V + dim W = 2d, then dim (VnW) = 
dim{V^nW^). 

Proof. Observe that (V + W)"^ = nW^. Since the subspaces V + W 
and {V + VF)" are symplectic orthogonal complements, the sum of their 
dimensions is 2d. Therefore, 

dim {V^ n VF^) = 2d- dim {V + W) 

= 2d-dimV - dimW + dim {V D W) 
= dim {VnW), 

which completes the proof. □ 

Theorem 2.15. Let G ^ P be a 2d- dimensional symplectic groupoid, and 
suppose T, C G is a Lagrangian submanifold. Then the restricted source map 
5|e : S — 7- P is a local diffeomorphism if and only if the restricted target map 
S^P is. 

Proof. For any fi £ T,, let V = kerT^S (i.e., the tangent space to the source 
fiber at //) and W = T^jT,. The source and target fibers are symplectic 
orthogonal complements, so = kerT^/S, and because the fibers have 
equal dimension, it follows that dimF = dimy^ = d. Next, because S 



10 



J. C. MARRERO, D. MARTIN DE DIEGO, AND A. STERN 



is a Lagrangian submanifold, this implies that = W, so dimVF = d. 
Therefore, diml/ + dimM^ = 2d, so these subspaces of T^G satisfy the 
conditions of lLemma 2.14[ 

Finally, ajs is a local diffeomorphism at fi if and only if the kernel 

ker (5|s) = ker T^S n T^S = VnW 

is trivial. Likewise, f3\-£ is a local diffeomorphism at fi if and only if the 
kernel _ _ 

kerr^(/3|s) = kerr^/3 n T^S = V^nW^ 



is trivial. But Lemma 2.14 implies that these kernels have equal dimension, 



so one is trivial if and only if the other is. □ 
Remark 2.16. Although we are only concerned with the conditions for 5|s 



and /3|s to be local diffeomorphisms, the proof of Theorem 2.15 also applies 
to any other property that can be described in terms of the dimension of the 
kernels of their tangent maps. (For example: 5|s has constant rank if and 
only if does.) 



Using Theorem 2.15 we deduce that if a|x; is a local diffeomorphism, then 
so is and hence the discrete Lagrangian flow map (5|x;)~"'^o/3|s is a local 
automorphism on S. Reversing the order of composition, it also follows that 
the discrete Hamiltonian flow map o (5|e)~ is a local automorphism on 
P — and moreover, it is a local Poisson automorphism. To see this, consider 
the Poisson map (5, /3) : G — )■ P~ x P, fj. and observe that 

the image of S is precisely the graph of /3|2 o (Sis)" in P~ x P- However, 
since S is Lagrangian, its image under the Poisson map (5, f3) is coisotropic; 
thus, it follows from a result of Weinstein [21] that /3|s o (5|s)~^ is a (local) 
Poisson automorphism on P, since its graph is coisotropic in P~ x P. This 
argument is essentially due to Ge |5j , who further showed that this map from 
(local) Lagrangian bisections of a symplectic groupoid to (local) Poisson 
automorphisms on the base is a group homomorphism. (The group structure 



that allows for composition of bisections is discussed in Appendix A 



3. Discrete constrained Lagrangian mechanics 

3.1. Generating Lagrangian submanifolds of T*G. In this section, we 
will be concerned with generating a Lagrangian submanifold C T*G 
of the cotangent groupoid, associated to a function L called the discrete 
Lagrangian. This is a particular example of the formalism introduced in the 



previous section. As shown in Example 2.13, one way to do this is by defining 



a discrete Lagrangian L: G — ?■ M and taking T,i = dL{G); this is the case 
considered in previous work on discrete Lagrangian mechanics, including 
Weinstein [22] and Marrero et al. [15] . 

Setting aside the groupoid structure for the moment, this approach exploits 
the fact that for any manifold M and function L: M — )• M, the submanifold 
dL(M) C T*M is Lagrangian. However, there is a more general construction 
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due to Tulczyjew (and stated, for example, in [20]), which we will use to 
generalize the earlier approach to discrete mechanics. 

Theorem 3.1 (Tulczyjew). Let M be a smooth manifold, N C M a sub- 
manifold, and L: — ?• R. Then 

Si = G T*M I TTMip) e N and {p, v) = (dL, v) 

for all V G TN C TM such that tm{v) = ttm{p)} 
is a Lagrangian submanifold ofT*M. 

(Here, ttm- T*M — )• M and tm- TM — )• M denote the cotangent and 
tangent bundle projections, respectively.) In the special case = M, this 
gives the familiar Lagrangian submanifold = dL(Af) C T*M. 

Turning back to the groupoid formulation, this allows one to define a 
discrete Lagrangian on some constraint submanifold N G G, rather than 
necessarily on all of G. This realization motivates the following definition. 

Definition 3.2. A discrete constrained Lagrangian system consists of a 
triple {G,N,L), where G ^ Q is a Lie groupoid, A^ C G is a submanifold, 
and L: A^ — )• M is a function called the discrete Lagrangian. 



It follows immediately from Theorem 3.1 that a discrete constrained 
Lagrangian system generates a Lagrangian submanifold 'Sl G T*G oi the 
cotangent groupoid T*G ^ A*G. The relationship among these spaces is 
shown in the following diagram: 

Sl c >T*G^^A*G 



N 



^G 



Definition 3.3. Let (G, A^, L) be a discrete constrained Lagrangian system. 
Then define the discrete Legendre transformations ¥^L: — )• A*G to be 
the restricted source and target maps, 

F+L 



W-L = 5|s, 



/5h 



Corollary 3.4. The discrete Legendre transformation ¥ L is a local diffeo- 
morphism if and only if¥^L is. 



Proof. Direct application of Theorem 2.15 



□ 



Definition 3.5. A discrete constrained Lagrangian system (G,N,L) is said 
to be regular if F^L are local diffeomorphisms, and hyperregular if they are 
global diffeomorphisms. 



Now, given a trajectory ^ui, . . . , /i„ G T*G, it follows from Definition 2.12 



that this satisfies the discrete Lagrangian dynamics of when fii, 
Si and 

F+L (;Ufc) = F-L (/.fc+i) , A; = l,...,n-1. 
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Hence, if {G, N, L) is regular, then the discrete Lagrangian flow map is 
given by the local automorphism (F^L) o F^L on S^, while the discrete 
Hamiltonian flow map is given by the local Poisson automorphism F^L o 
(F~L) ^ on A*G. In the special case N = G, T^l = dL(G), this is in 
agreement with the formulation of Weinstein |22j and Marrero et al. |i3j, as 



discussed in Example 2.13 



Remark 3.6. It should be emphasized that this formalism does not require 
to be a subgroupoid of G, but only a submanifold. If C G is indeed 
a subgroupoid, then one can simply reduce to the unconstrained dynamics 
on ^ Q, which are given by the Lagrangian submanifold d-L(A^) C 
T*N ^ A*N. Therefore, in the subgroupoid case, N can be thought of as 
specifying discrete holonomic constraints on G. By contrast, when is not 
a subgroupoid, it can be thought of as specifying discrete nonholonomic 
constraints on G. 

This is consistent with the definition of holonomic and nonholonomic 
constraints for continuous Lagrangian systems. For example, a constraint 
distribution A C TQ is integrable precisely when A is closed under the 
Jacobi-Lie bracket on TQ, i.e., when A is a subalgebroid of the tangent Lie 
algebroid TQ. Just as continuous holonomic constraints correspond to Lie 
subalgebroids, it is natural to think of discrete holonomic constraints as 
corresponding to Lie subgroupoids. 

3.2. AfRne bundle structure of S/,. For any discrete constrained La- 
grangian system {G,N,L), the Lagrangian submanifold C T*G is also a 
bundle over A^. More precisely, taking the projection to be the restriction 
of the cotangent bundle projection itg\t:l- A^, we obtain an affine 

bundle whose associated vector bundle is iy*N, the conormal bundle of A^ in 
G. To see this, note that for any G and AG i/*A^ at the same basepoint 
TTGifJ-) = ttg{A-) = g, we have 

{fi + A,v) = = {dL{g),v) 

for all V G TgN, and thus /i + A G S^,. (For more details, see Libermann and 
Marie [TU] and references therein.) 

Therefore, is isomorphic to i'*N, but generally not in any canonical 
way. To choose a particular isomorphism, if one so desires, it suffices to 
specify a distinguished section a : N ^ T^l, and then 

= Wig) + A I A G iy*N, g = vtgIA)} ^ u*N. 

In particular, suppose that L is an extension of L to a neighborhood of A^ in 
G, so that L = L\n- Then this defines the distinguished section a = dL\i\i, 
and hence 

Si = {dL{g) + A I A G u*N, g = ^g(A)}. 

Applying the definitions of a and /3, it follows that the discrete Legendre 
transformations can be written as 

(F-L(dL(5) + A),X(a(5))> = {dL{g) + A,l (g)) = l[L]ig) + (A, ^(5)), 
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and similarly, 

{¥+L{dL{g) + A),X{P{g))) = ^[L]{g) + (A, 1({g)), 

for all sections X £T (AG). 

Now, given a sequence G T*G, let gk = t^g it^-k) and take 

= fJ-k — dL (gk) for = 1, . . . , n. Therefore, this is a solution of the 
discrete Lagrangian dynamics when Ai, . . . , A„ ^ v*N and 

(5) ^[L] + (A,., ^ (5fc)> = (5fc+i) + (Afc+i, ^ (5fe+i)), 

for all sections X {AG) and A; = l,...,n — 1. In the special case N = G, 
observe that i'*N is simply the zero section of T*G, which is isomorphic to 
G itself; hence, the dynamics reduce to the unconstrained discrete Euler- 
Lagrange equations of Example 2.13[ 

3.3. Lagrange multipliers. Suppose that the constraint submanifold A'^ C 
G is defined by 

N = {gGG\ r{g) = 0, a G ^} , 

where {^'"'^a&A ^ family of real functions defined in a neighborhood of N 
and A is an index set. It follows that ^^j^^ is a basis of sections of 

the conormal bundle i'*N . Hence, a section A of the conormal bundle can 
be written 

A = Aad^'^l^, 

where the real functions Aa on N are called Lagrange multipliers. (Here, we 
have used the Einstein summation notation to indicate that we are summing 
over a ^ A.) In fact, since 0"|Ar = 0, we can deduce that 

A = d(Aa0l7v). 

Now, as before, suppose that the discrete Lagrangian L: — t- M is the 
restriction to A^ of a real function L on G. Therefore, an element /x G 
with g = TTcif^), can be written as 

M = dL{g) + d (Xa^) (g) = d(L + A^c^") (g) G E^. 

In this sense, Til can be seen as the space consisting of elements g £ N, 
together with the Lagrange multipliers A constraining g to A^. 
Thus, if {g, X) £ T,L, we have 

{¥-L {g, A) , X (aig))) = (d(L + Xa^) (g), 1(g)) = 1[L + Xa^] (5), 
and likewise 

(F+L (9, A) , X if3{g))) = ^ [L + A,^'^] (g), 

for all sections X £ T (AG). Consequently, let (gi, Ai) , . . . , {gn, A„) be a 
sequence of groupoid elements and Lagrange multipliers. Then this is a 
solution of the discrete Lagrangian dynamics when gi, . . . , gn £ N, i.e., 

0" idk) = for all a G ^, fc = 1, . . . , n. 
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and when 

^ [L + (Afc), r] {gu) = 1[L + (Afc+i), <p''] (gk+i) , k = l,...,n-l, 
for all sections X £T (AG). 

4. MORPHISMS, REDUCTION, AND NOETHER SYMMETRIES 

4.1. Groupoid morphisms and reduced dynamics. In order to study 
reduction of discrete constrained Lagrangian systems, we first recall the 
definition of a morphism of Lie groupoids. After this, we introduce the 
slightly more specialized notion of a morphism for discrete constrained 
Lagrangian systems, which preserves not only the groupoid structure, but 
also that of the constraint submanifolds and the Lagrangian functions. 

Definition 4.1. Given two Lie groupoids, G ^ Q and G' ^ Q' , a smooth 
map <I> : G — )• G' is a morphism of Lie groupoids if, for every composable pair 
{g,h) £ G2, it satisfies {<S>{g),^{h)) £ G'^ and $ (gh) = <^>{g)^{h). 

Definition 4.2. Given two discrete constrained Lagrangian systems, denoted 
{G, N, L) and {G' , N' , L'), a smooth map $ : G — )• G' is a morphism of discrete 
constrained Lagrangian systems if it is a morphism of Lie groupoids, and 
additionally, it satisfies N = (iV') (i.e., g e N 4^ $(5) e A^') and 
L = L' o 

Example 4.3. Suppose that {G,N,L) and {G' , N' , L') are discrete con- 
strained Lagrangian systems, with constraint submanifolds defined by 

N = {gGG\ <l>''ig) = 0, a e A} , N' = {g' e G' \ ^ (5) = 0, a G ^} . 

If G — )• G' is a morphism of Lie groupoids satisfying = 0'" o then 
clearly 4>°-{g) = 4^ 4>"^ {^{9)) = 0, so g e N 4^ $(5) G N'. If, furthermore, 
we have L = L' o $, then this implies that $ is a morphism of the discrete 
constrained Lagrangian systems. 

A morphism $ : G — )• G' of Lie groupoids induces a smooth map ^q: Q ^ 
Q' on the base, which satisfies 

a' o $ = $0 o a, /3' o $ = c|)q o /3, <I) o e = e' o $0, <I> o i = o $0 . 

Moreover, ^> induces a morphism A^ : AG — )• AG' of the corresponding Lie 
algebroids, and 

'^Hv){g) = r4(3) {A^{v)) = TcD {Tlgiv)) = TcD (^(5)) , 

AHw){g) = -T (r$(g) o i') {A<^{w)) = T<^ {-T {rg o i) (w)) = (li{g)) , 

for any g £ G and any v G 74^(g)G, w G A^^^g-^G. Consequently, two sections 
A G r (AG) and A' G L {AG') satisfy A$ o A = A' o $0 if and only if 
= A' o <!>, or equivalently, r$ o = A' o $. (See, e.g., Mackenzie 
p. 125].) That is, two sections are "^$-related" if and only if their 
corresponding left-invariant (and right- invariant) vector fields are ^-related. 
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tA*G' 




Figure 1. A diagrammatic depiction of Theorem 4.5, show- 
ing the maps and relations involved in the reduction of discrete 
constrained Lagrangian systems. 

It is now possible to define the dual notions of these relations on the cotangent 
groupoid and its base. 

Definition 4.4. Two covectors fi £T*G and fi' £ T^^^^G' are ^* -related if 
{fi,v) = {fi',T<^{v)) for all v G TgG. Similarly, p G A*gG and p' G Al^^^^G' 
are A* <P -related if {p,v) = {p' ,A^{v)) for all v G AqG. 

We now have the necessary equipment to state and prove the main theorem 
on the reduction of discrete constrained Lagrangian systems, which states 
how the structure of the cotangent groupoid and Lagrangian submanifold 
transform under morphisms. As an immediate corollary, we see that the 
Lagrangian dynamics of (G, A^, L) reduce to those of (C, A', L'), with respect 
to $*-relatedness. (The content of this theorem is shown diagrammatically 
in 



Figure 1 



Theorem 4.5. Let G ^ G' he a morphism of the discrete constrained 
Lagrangian systems {G,N,L) and {G' , N' , L'), and suppose // G T*G and 
fi' G T*G' are ^* -related. Then, the following are true: 

(i) if n' G Sl', then ii G S/,; 

(ii) the sources a{^) G A*G and a' {n') G A*G' are A* ^ -related; 

(iii) the targets e A*G and P' {fJ) G A*G' are A* ^-related. 

Proof. To prove (1), we begin by denoting g = TrcifJ') and g' = ttq' {fi'); since 
fx and fi' are $*-related, we have ^{g) = g' . If ^' €T,l/, then ^{g) = g' e A', 
so 5 G (A') = A. Furthermore, for all v G TgN, 



[dL',T^v)) 

= (cD* (dL') , = (d {L' o<^),v) = (dL, v) , 
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and thus fi gT,l. 

To prove (2), suppose that v G Aa(^g)G. Then 



so a{fi) and a' {fj,') are ^*<I>-related. Likewise, to prove (3), let v G Af^f^g-jG. 
Then 



(/3(/x),^> = (a., V> = (/.',TcI>(V)> = {fs',A^v)) = (/?' {fi'),A^v)), 
and thus and /3' {fi') are ^*<^-related. □ 

Corollary 4.6. Let $ : G — )• G' be a morphism of the discrete constrained 
Lagrangian systems (G, A^, L) and {G' , N' , L'). If ii'^, . . . , G T*G' satisfies 
the discrete Lagrangian dynamics of {G' , N' , L'), then any ^* -related sequence 
fii, . . . , fin G T*G satisfies the discrete Lagrangian dynamics of (G, N, L). 



Proof. By Theorem 4.5, /i^ S Ti^i impHes /i^ G S^, for A: = 1, . . . ,n. Further- 
more, for any v S Ap(^g^^G = ^^.(g^^^jG, we have 

(p{lik),v) = {0 {^x'^),A^v)) = {a' (a^Ui) ,A^{v)) = {a{fik+i),v), 

so /3 (/ifc) = 5 (/xfc+i) for /c = 1, . . . , n — 1. Thus, //i, . . . , fin satisfies the 
discrete Lagrangian dynamics of {G,N,L). □ 

Remark 4.7. If $ : G — )■ G' is a submersion, then the converse is true as well. 
Since every v' G T^(^g-^N' can be written as v' = T^{v) for v £ TgN, part (1) 
of Theorem 4.5 can be strengthened to /x G 44> /x' G S^/. Furthermore, 
A^ is a fiberwise surjection if (and only if) $ is a submersion (Mackenzie 
im Proposition 3.5.15]). Therefore, with this additional assumption, Corol- 
|lary 4.6 can be strengthened to say that the sequence fii, . . . , fin G T*G 
satisfies the discrete Lagrangian dynamics if and only if the -related 
sequence fi[, . . . , fi'^ £ T*G' does. 

4.2. Noether symmetries and constants of the motion. Next, we 
extend the notion of discrete Noether symmetries and constants of the 
motion, as defined for unconstrained systems on Lie groupoids in Marrero 
et al. 113], to apply to discrete constrained Lagrangian systems. 



Definition 4.8. A section X G F (AG) is said to be a Noether symmetry 
of the discrete constrained Lagrangian system (G, A^, L) if there exists a 
function / G G°°((5) such that 

{aifi),X {a{g))) + f {a{g)) = {P{fi), X i(3{g))) + / (Pig)) 

for all ;U G Sl, where we denote g = -Kcifi)- 

For each Noether symmetry of a discrete constrained Lagrangian system, 
there is a corresponding constant of the motion, which is preserved by the 
Lagrangian dynamics; this is the discrete version of Noether's theorem. 
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Theorem 4.9. If X £ T (AG) is a Noether symmetry of the discrete con- 
strained Lagrangian system (G, N, L), then the function Fx : S_l — )• M defined 
by 

= (5(/x), X {a{g))) + f {a{g)) = 0{fi),X {(3{g))) + / {P{g)) , 

where g = ■kq{^), is a constant of the motion. That is, if . . . , fj,n £ T*G 
satisfies the discrete Lagrangian dynamics, then Fx (fJ-k) = Fx (/^fc+i) for 
k = 1, . . . ,n — 1. 

Proof. If /xi, . . . , /i„ G T*G satisfies the discrete Lagrangian dynamics, then 
by definition, fi^ G for k = l,...,n and (3 (fik) = OiifJ-k+i) for k = 

1. . . . , n — 1. Therefore, 

Fx (//fc) = 0{f^k) , X (/3 (5fe))> + / (/? {gk)) 

= (5 {jxk+i) , X (a {gk+i))) + f{a (gk+i)) 
= Fx (pk+i) , 

for k = 1, . . . ,n — 1, which completes the proof. □ 

5. Variational principles 

Discrete constrained Lagrangian systems are a generaHzation of variational 
integrators. However, up to this point, we have only treated the discrete 
Lagrangian L: A^cG— )-Masa generating function, and have not yet made 
any mention of its variational interpretation. In this section, we show that, 
under certain regularity conditions, the Lagrangian dynamics of (G, N, L) 
are precisely the critical points of a discrete constrained action sum. 

5.1. Admissible trajectories and variations. Given a fixed element g € 
G, a trajectory gi, . . . ,gn £ G is said to be admissible if gk £ N for k = 
1, . . . , n, and if the elements are composable with product gi ■ ■ ■ gn = 9- Let 
^gN denote the set of admissible trajectories of length n. Without loss of 
generality, we can take n = 2, so that the space of admissible trajectories is 

Cg,N = {(51,52) G G2 n (iV X iV) I gig2 = g} = {mlG^niNxN))''^ io)- 

Finally, assume that G2 n {N x A^) is a submanifold of N x N, and that 
the map m\G2n{NxN) '■ G2 H (A^ x A^) — > G has constant rank in an open 
neighborhood of C|^. This is sufficient to ensure (by the subimmersion 
theorem, cf. Abraham et al. pj) that is itself a submanifold of G2 n 
{N X A^). 

Now, a "variation" of an admissible trajectory (51,52) ^ C^at is an element 
of the tangent space 

'^{gi,g2)^lN = {{Vg^.Vg^) £Tg^NX Tg.2N \ Tg^^ {Vg,) = Tg^a [Vg^) 

and T^g^^g.2)m{Vg^,Vg2) = O}. 
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In fact, letting q = (3 (gi) = a (52), this tangent space is isomorphic to a 
vector subspace {AqG)^j^ ^^-^ C AqG, defined by 

(AG)(;v,,i,32) = e AG I V (51) G Tg,N and it (52) G Tg,N} , 

where we recall that V {gi) = T^^^g-jig-^^v) and it {g2) = —T^{q) (fg^ o i) (v). 
The following theorem establishes this isomorphism explicitly. 

Theorem 5.1. The linear map ip{N,gug2)- i^Q^)iN,gi,g2) ~^ '^{9i,92)Cg,N> 
defined by v ^ (V {gi) , —it (32)), is an isomorphism. 

Proof. The proof proceeds in three stages: first, showing that f(N,gi,g2) 
well-defined; second, that it is injective; and finally, that it is surjective. 

¥'(Af ,91,92) well-defined: Given any v G {AqG)^j^ we wish to show 
that (/'(Af,gi,g2)(^) ^ '^(9i,92)Cg,N- % definition of the space {AqG){N,gug2)' 
clearly have (V (gi) , -it (52)) G TgjiV x Tg^N. Moreover, 

r/3(V(5i)) = p(t;) = Ta(-lt(52)), 

so all that remains is to show that Tm (V (171) , —it (52)) = 0. Using the 
results from Appendix A let us suppose that Bi and B2 are local bisections 
of G such that gi G Bi and (72 G i?2, so (-Bi)^ (g) = 51 and (-62)0, (q) = 52- 
Then, applying Theorem A.l and the inversion identity (|3j), we obtain 

Tm (V (gi) , -it (52)) = Tm (T^^, (i;), T (r^, o i) (t;)) 

= T {rB2 ° ^Bi ){v) + T {Ib^ o rB2 ° i) {v) 

- T o o e o /3) (v) 
= T {rB2 o iBi ) {v) - T {(bi o rB2 ) {v) 
which vanishes since ° = C-Bi ° fB2- Therefore, we have shown that 
'P(N,gi,g2)iv) G r(g^,g2)C2^, as claimed. 

^(N,gi.g2) '^s injective: This follows immediately from the fact that the 
linear map T^{q)igi ■ AqG — )• Tg-^G, w 1— V (gi), is injective. 

'f{N,gi,g2) is surjective: Suppose that (^91,^92) G T(^g-^ g^^Cg j^. Then we 
can take some curve in Cg j^, 

c: {-6, 6)CR^ CIn, t ^ c{t) = (ci(t), C2(t)) , 

such that c(0) = (ci(0), C2(0)) = (51,52) and c(0) = ivg^,Vg.^). 
Now, since a (ci(t)) = a (gi) for all t, we can write 

ci{t) = gih{t), 

where h{0) = e{(3{gi)) and h{t) = g^^ci{t) G a^^ (Mdi))- Therefore, this 
defines a Lie algebroid element v = h{0) G AgG, and 

7;,, =ci(o) = r£,,(A(o)) = V(5i). 

On the other hand, since ci{t)c2{t) = gig2 for all t, it follows that C2{t) = 
h{t)~^g2. Thus, 

Vg, = C2(0) = T {rg, o i) (A(0)) = -It (52) . 
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Therefore, we have shown that v E (^<?^)(Af gi 92) that {vg^^Vg^) = 
V(Ar,gi,g2)(^)' f{N,gi,g2) surjective. This completes the proof. □ 

5.2. The discrete constrained action principle. Now that we have char- 
acterized the admissible trajectories of length two, and their variations, we 
introduce the discrete constrained action sum 



S{G,N,Ly -.CIj^^R, {gi,g2)^L{gi) + L{g2) 

-'9,N 



A critical point of S {G, N, L)"^ is a trajectory (51,52) £ such that 



dS{G,N,Lf{gi,g2)=0, i.e., 

= {dS (G, N, Lf , (V (51) , (52))) = (dL, V (5i)> - (dL, 1^ (52)) , 

for all f G where as before q = {3 (gi) = a (52)- 

To express these discrete constrained Euler-Lagrange equations in terms of 
arbitrary sections X gT (AG), as with the earlier results, suppose now that 
L is an extension to a neighborhood of C G, so that L = L|jv. Therefore, 
it follows from the above that 

(6) d(Loig^+Lorg,oi){e{q)) 



1 







= T:^g/9^ (d^ (51)) + T:^g) {Tg, O i) (dL (52)) G ( A^) J^,,, , 

where T*^^^£g^ and T*^^ [rg^ o i) denote the cotangent lift (i.e., the adjoint to 
the tangent map) of ig-^ and rg^ o i, respectively, and where (A^G)^^^^ C 
T*^g^G is the annihilator of (^qG)j^^^^ Now, by using the definition of 
(AqG)^]^ we can write its annihilator as 

(AG')(7V,gi,32) 

= Trgi {Tg, NnkevTg,a)nT {rg, oi)-^{Tg,Nn ker Tg, (3) 
= [Trg^ {Tg.Nn keiTg.a)]^ + [t (rg, o i)-^ {Tg,N n keiTgJ) 
= ^g)^9i {{Tg,Nn kerTg,af) + T:^^^ {rg, o i) ((r,,iV n kerTgjf 
Therefore, from ([6|, we obtain that 

T:^g^eg, (dL (51) + Ai) = -r;(^) (r,, o i) (dL (52) + A2) , 

where we have denoted Ai G {Tg-^^ N nkei Tg-^a)^ = (Tg-^^N)^ + (ker T^^a)'^ 
and A2 G (Tg^N D kei Tg^pf = (Tg^Nf + {ker Tg,(3f. Consequently, we can 
write 

Ai=Ai+Ai, A2 = A2 + A2, 
where Ai G {Tg^Nf = i^*g^N, M G (ker T^.a)", A2 G (r<;2A^)° = ^^^2^' and 
A2 G {keiTgjf. Finally, since Ai annihilates a- vertical tangent vectors, in 
particular it annihilates any left-invariant vector field evaluated at gi; likewise. 
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A2 annihilates any right-invariant vector field evaluated at g2- Therefore, it 
follows that for any X gT (AG), we have 

^[L] (51) + (Ai, ^ (gi)) = ^[L] (52) + (A2, 1 (52)), 

which is precisely the equation ^ that we derived earlier, in the case n = 2. 

To generalize this argument for any n > 2, one can proceed in a similar 
manner. Given a fixed g G G, suppose that the space of composable sequences 
of length n in N, 

= {{gi, ...,gn)GN \ {gk,gk+i) G G2 for A: = 1, . . . , n - 1} , 

is a submanifold of A^", and that the multiplication map 

m"|Arn : N2 ^ G, {gi, ■ ■ ■ , gn) ^ gi ■ ■ ■ gn, 

has constant rank in an open neighborhood of C^j^ = (m."|Ar^i) ^ (g). Then, 
as before, it follows that C^^y is a submanifold of N2. The discrete constrained 
action sum is then 

n 

S {G, N, Lf ■■ C^j, ^ M, {gi,...,gn)^Y.^ ^9k) , 

fc=i 

and one may prove that when {gi,...,gn) G is a critical point of 

S (G, A^, L)", the equations 

1([L] {gu) + (Afc, ^ {gk)) = 1[L] (gk+i) + {Ak+i,l (gk+i)), 

hold for all sections X gT {AG) and /c = l,...,n — 1, where A^ S '^^^ for 
k = 1, . . . , n. Again, this agrees precisely with the earlier equation ([5]). 

6. Examples 

6.1. Constrained mechanics and optimal control on Lie groups. Let 

G be a Lie group, whose identity element is denoted by e G G, and let q be 
the Lie algebra of G. Suppose that a submanifold A^ C G is given by the 
vanishing of the constraint functions (j)"", for a £ A, defined on a neighborhood 
of A^. Let L: A^ — )• M be a discrete Lagrangian, and let L be an extension of 
L to a neighborhood of A^, so that L = L\]\f. Then the discrete dynamics 
are given by the constraint equations 

idk) = for ah a G ^4, k = 1,... ,n, 

together with the equations 

(7) t[L + {Xk)a r] {9k) = t [£ + (Afc+i), r] {gk+i) , k = l,...,n-l, 
for every C G 0- 

Now, in addition to the usual definition ^ {g) = T£g{^), we can use the 
tangent inversion identity ^ to see that Ti{S,) = T {e o f]) (^) = — ^; 
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therefore, by the chain rule, we simply have ^ (g) = —T (r^ o i) = Trg{^). 
Hence, can be rewritten as 

l^d (l + (Afc), r) (e) = 4^,d {l + (Afc+i), r) (e), A; = 1, . . . , n - 1. 

If we define /Jk = ^gfcdL(e) G g* and = r*^d(/)"(e) e g* for each k, then 
this is equivalent to 

Mfc+i + {h+i)a ^k+1 = Ad*^. {f^k + {^k)a ^fc) > A: = 1, . . . , n - 1, 

where Ad* : G x g* — )• g* is the coadjoint action of G on g*. These equations 
will be called the discrete constrained Lie-Poisson equations for this system. 

We now show how this framework can be used to discretize a general 
family of optimal control problems on Lie groups. Again, let G be a Lie 
group with Lie algebra g, and define some control submanifold G C g. A 
curve ^ : [0, 1] — )■ G is said to be a control curve for a trajectory g: [0, 1] — )• G 
if it satisfies g{t) = Tigi^^-^ {^{t)) for all t; when G is a matrix Lie group, we 
can simply write this as g{t) = g(t)^(t). The continuous optimal control 
problem is defined as follows (cf. Koon and Marsden [9]): 

Given an initial configuration go £ G, a final configuration 
gi £ G, and a function [: g — >• M, find the control curve 
^: [0, 1] — > G for a path g: [0, 1] G, satisfying g{0) = go 
and g{l) = gi, such that dt is minimized. 

In other words, we wish to find the most efficient (lowest-cost) way to "steer" 
the system on a path from gQ to gi by applying controls in G. (In optimal 
control problems, it is generally assumed that the system is controllable, i.e., 
that there exists a control curve for every pair of initial and final configu- 
rations. As in Koon and Marsden [9j, we omit discussion of controllability, 
since we are primarily interested in the necessary variational conditions for 
optimality.) 

To derive the associated discrete system, it is usually necessary to introduce 
an analytic local diffeomorphism r : g — ?■ G, which maps a neighborhood of 
G g to a neighborhood of e € G, where e denotes the identity element of G. 
As a consequence, it is possible to deduce that T(^)r (— ^ = e. This map can 
then be used to relate the continuous and discrete settings (cf. Bou-Rabee 
and Marsden [2]). There are several choices, in the literature, for the map 
r. For instance, one may take r to be the exponential map, defined by 
the time-1 fiow exp(^) = 7(1), where 7: [0, 1] — t- G is the integral curve of 
the vector field G X(G) with initial condition 7(0) = e. Alternatively, 
for quadratic matrix Lie groups (e.g., S0(3), SE(2), SE(3)), it is typical to 
use the Cayley map cay: g — )• G, defined by cay(^) = (/ — ^/2)^^ (I + C/2), 
where I = e denotes the identity matrix. 

Now, suppose that we start with a continuous optimal control problem, 
specified by a Lagrangian [: g — )• M and a control submanifold G determined 
by the vanishing of the independent constraints : g — )• M, a G j4. Then, 
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given a choice of the map r : g — )• G, we may construct the discrete Lagrangian 
L : G — >• M and constraint functions (j)"" : G ^ M, a G A, hy 



(r-io£,J*(d[(6) + (A,),dM/'^(a.)) 

= {r~' ° rg^^,)* (d[ (Cfc+i) + {Xk+i)a fe+i)) , k=l,...,n-l, 



where = {gk) /h. In other words, denoting dr^ : g — t- g to be the 
left-triviahzed tangent map of r, defined by T^t = Telr(£) ° dr^, the previous 
equation can be rewritten as 



(dT,/J*(d[(6) + (Afc)„d^'^ (6)) 

= J* (dl (6+1) + fe+i)) , = 1, . . . , n - 1, 



6.2. A discrete plate-ball system. Consider, now, the optimal control 
problem for the classical example of a homogeneous ball, rolling (without 
slipping) on a rotating plate. The configuration space of this system is 
Q = xS0(3), parameterized by g = (x, y, g). Let (w^., w^, ojz) G M'^ = 50 (3) 
be the angular velocity vector of the ball with respect to an inertial frame, 
and let r be its radius. Suppose that the plane rotates with constant angular 
velocity about the axis perpendicular to the plane through the origin. 
Then, we may formulate the following optimal control problem, which is 
called the plate-hall problem: 




where h is the time step size of the discretization. 
Therefore, we can finally rewrite Equation ([T]) f 



as 



since we have dr^ = T^{^)T o Teir(^)- 



Given initial and final configurations qo,qi G M? x S0(3), 
find the optimal control curves {x{t),y{t)) G M^, t G [0, 1], 
that steer the system from qq to qi, such that the cost 
function 




is minimized, subject to the (nonholonomic) rolling-without- 
slipping constraints 



y + rujx = 



X — rujy = —Vty 



UJz = C. 



where c is a constant. 



The continuous equations of motion for this problem are carefully studied 
in Koon and Marsden [9|, Iglesias et al. [7j; we now study the geometric 
discretization of this system and its dynamics. 
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The Lie groupoid that arises for the plate-bah problem is ] 



'xS0(3) 



M?, where the structure functions are given by 

a{xo,yo,xi,yi,gi) = (xo,2/o) , 
f3 {xo,yo,xi,yi,gi) = {xi,yi) , 
m{{xo,yo,xi,yi,gi) , (xi, yi, X2, y2, 52)) = ixo,yo, X2,y2, 9192) , 
i {xo,yo,xi,yi,gi) = (xi, yi, xq, yo, ^f^) ' 
e{x,y) = {x,y,x,y,I) . 

Now, denote by {Ei, E2, E^} the standard basis of the Lie algebra so(3), 



E^ 




Eo 




E-, 




so (3) (sometimes called 






-UJz 


UJz 





-UJy 





This defines the Lie algebra isomorphism ^: W 
the "hat map" [I]), which takes 

UJ = {ujr,, UJy,UJz) '-^ U) = UJxEi + UJxE2 + UJxEs 



As a discretization procedure for the Lie group S0(3), it is typical to use the 
matrix logarithm, denoted log: S0(3) — )• so(3), which is the (local) inverse 
of the matrix exponential map exp: so(3) — t- S0(3). For simplicity, we may 
use the approximation exp (/i^) / + /i^ to obtain 

logy 9- 1 



h 



Taking tD = ^, we therefore deduce that 

1 



UJx 



UJv 



UJ, 



trim 
tr {CE2 
tr(Ci?3 



h 



1 

'2h 
1 

'2h 
1 

'2h 



tr {gEi 
tr {gE2 
tr (y^3 



since tr {lEi) = tr {Ei) = for i = 1, 2, 3. 

Finally, we derive the following discretization of the plate-ball system. Let 
N be the submanifold of SO (3) determined by the vanishing of 

the constraint functions, 



(j)^ (xo,yo,xi,yi,yi) = h 
(8) 0^ (xo,yo,xi,yi,yi) = /i 
(j)^ (xo,yo,xi,yi,yi) = h 



yi - yo 

h 

xi - Xq 



h 



-^tr{g,E,)-n^-^±^ 
+ ^tr{g,E2) + n^^ 



1 



c+^tr(yii?3) 
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and take the discrete Lagrangian to be 



L{xo,yo,xi,yi,gi) 



h 



+ 



y\ - yo 

h 



Therefore, applying the results of Section 3.3, the discrete constrained La- 
grangian dynamics are given by 



= 1, . . . , n — 1, 



+ {\k)a<P"'] {xk-i,yk-i,xk,yk,gk) 
= ^[l+ {\k+i)a(l^"'] {xk,yk,Xk+i,yk+i,9k+i] 

together with the vanishing of the constraint functions ([s]), where X is an 
arbitrary section of the vector bundle TM^ x so (3) — )■ M^. A basis of sections 
of this vector bundle is given by 

Hence, the discrete dynamics are given by the following system of equations: 



Xk+i — '2xk + Xk-i ^ (Afc+i)2 



/l2 



h 








- 2yfc + , (Afe+i)2 - (Afc)2 ^{^k+i)i + i^k)i 



+ 









/l2 ' h 

-r (Afe)i tr {gkEf) + r {Xk+i)i tr {Eig^+iEi) 
+ r (Afc)2 tr {gkEiE2) - r {Xk+i)2 tr (£^i5fc+i^2) 

- (Afc)3 tr {g^EiEs) + (A^+Og tr (^15^+1^3) , 
-r {Xk)i tr igkE2Ei) + r {Xk+i)i tr (£^25fc+i-Ei) 

+ r {Xk)2 tr {gkE2) - r {Xk+i)2 tr {E2gk+iE2) 

- (Afc)3 tr {gkE2E-i) + {Xk+i)^ tr {E2gk+iEz) , 
-r (Afc)i tr (gkEsEi) + r (Afc+i)i tr (E's^fc+i-E'i) 

+ r {Xk)2 tr {gkE-iE2) - r {Xk+i)2 tr (-^35/0+1^2) 

- (Afc)3 tr (c/fc^l) + (Afc+i)3 tr (£^35^+1^3) , 
yk+i - yk 



h 



W TP \ f^Xk+l+Xk 

— tr [gk+iEi) - n , 



— Xk r , ^ , ^yk+^ + Vk 
h 2h + ^ 2 

c + ^ tr (gk+iEs) . 



These are eight equations (corresponding to the five basis elements and 
three constraints), and each step of the dynamics requires solving for eight 
unknowns (the five degrees of freedom x, y, iVx, tOy, lOz, plus the three 
Lagrange multipliers). 
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6.3. Time-dependent constrained mechanics. To treat time-dependent 
discrete mechanics on a Lie groupoid G ^ Q, we construct the new Lie 
groupoid Gk = MxMxG^Mx(3, with structure functions defined by 

"R ito,ti,gi) = (to, a (51)) , 
^R {to,ti,gi) = (ti,/3(gi)) , 
"T-R ((^0,^1, 51) , (^1,^2, 92)) = (to, ^2, 5152) , 
«R ito,ti,gi) = (ti,to,5r^) ' 
eR {t,q) = {t,t,e{q)) . 

The associated Lie algebroid is naturally identified with TM x AG — t- M x Q. 
Hence, the sections of this Lie algebroid are spanned by O) and elements 
having the form (0, X), where X €T (AG). Therefore, we can express the 
left- and right-invariant vector fields as follows: 



,0) {to,ti,gi) - (-i|t=i^,0ti,0gi) -: -m\(to,tu9i)' 

(0,X){to,h,gi) = (0t„,0t,,^(5i)) =: ^(to,ti,5i) 
(0,X) (to,ti,5i) = (Oto,Ot„^ (gi)) =:1 {to,h,gi) 



Now, suppose we are given a constraint submanifold Nr C Gm. and a 
discrete Lagrangian : — t- M. Then, following the theory presented in 



Section 3, the discrete dynamics correspond to the Lagrangian submanifold 



T,!^^ of the cotangent groupoid T*G]r ^ A*Gm.- More explicitly, take an 
arbitrary extension Lr, and suppose that is given by the vanishing of 
constraint functions (p^, a £ A. Then the discrete dynamical equations are 

4>R{tk,tk+i, gk+i) = 0, a e A, 

and 

^[Lr + {^k)a'pR] {tk~l,tk,gk) + mi-^^ + i^k+l)a4>R] (4 , *A:+1 , ) = 0, 

1([Lr + {Xk)a^i\ {tk-i,tk,gk) -J^[Lr + (Afc+ l)a<AR| {tk,tk+l,gk+l) =0, 
for k = 1, . . . ,n — 1. 

Example 6.1. Let us construct an integrator for the constrained optimal 
control problem on a Lie group, as in Section 6.1, but incorporating adaptive 
time-stepping, rather than a fixed time step h. Define the time-dependent 
discrete Lagrangian 

-^^R {tk-l,tk, gk) = {tk — tk~l) I ( 

\tk — tk-1 
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along with the constraint functions 

(Pm {tk-i,tk,9k) = {tk - tk-i) (7^4^) ' « ^ ^• 

\i'k — t-k-l J 

Taking ^fc = r^-*^ {gk) /h^, where hk = tk — tk-i, the discrete dynamics are 
thus given by the equations 

= (Cfe+i) , 

0={dT^^\Sidli^k) + {Xk)ad^^m 

= [ (Cfc) - (dt (6) , - I (6+1) + {dl (a+i) , 6+1) • 

The last condition, corresponding to the ^ terms, can be interpreted as 
conservation of energy along the discrete evolution, as with the symplectic- 
energy- momentum preserving methods of Kane et al. [8j. 

Example 6.2. The previous example allowed for variable time steps, with 
constraints placed only on the Lie groupoid G. On the other hand, we may 
consider time-dependent discrete Lagrangian mechanics which are uncon- 
strained on G, but with fixed time step size tk — tk-i = h for k = 1, . . . ,n. 
As before, consider the groupoid Gjr =| M x Q, and define the constraint 
submanifold 

= {{to,ti,gi) £ Gr\ ti = to + h} , 

for some constant h. This corresponds to the vanishing of the single constraint 
function 

4>R {tQ,ti,gi) =ti-to-h. 

Then, given an extended discrete Lagrangian Lk : — )• M, the equations of 
motion are (together with the constraint equation) 

^ [Lr + Xk4>R] {tk~i,tk, gk) + §1 \Lr + \k+i(t)m\ {tk, tk+i,gk+i) = 0, 

^[-^M + Xk4>R] {tk-i,tk, gk) - X[Lr + Xk+i(j)M] {tk,tk+i, gk+i) = 0, 
for k = l,...,n — 1. However, this can be simplified greatly, since we observe 
that |, [^r] = 1, I [M = -1, and ^ [c))r\ = 1 [0m] = 0. Therefore, the 
last two equations become 

^[-^^r] {tk-i,tk, gk) + Afc + ^[-^^r] {tk,tk+i, gk+i) — A/fc+i = 0, 

^ [Lr] {tk-i,tk, gk) — ^ [Lr] (tfc, tfc+i, Sfc+i) = 0. 

Finally, observe that these two equations are completely decoupled, so we 
can in fact eliminate the first equation. Therefore, we obtain 

^[-^m] {tk-i,tk, gk) — ^[-^^r] {tk,tk+i, gk+i) = 0, 

with tk+i = tk + h, for A; = 1, . . . , n — 1. 

This has precisely the same form as the discrete Euler-Lagrange equations 
in the time-independent case; as in the continuous theory of Lagrangian 
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mechanics, no "^"-type terms arise from introducing the time dependency. 
In the special case where Lr depends only on the time step size tk+i — tk, this 
is equivalent to defining the time-independent discrete Lagrangian L : G — )• M, 
and one recovers precisely the usual, unconstrained discrete Euler-Lagrange 
equations. 

Finally, we remark that an extension of this setup can be used for more 
sophisticated step size control, by taking the constraint function to be (for 
example) , 

4>R ito,ti,9i) =ti-to-h (gi) , 
where h: G — )■ M is some step size function. In this case, ^ [(j)^] = [h] and 
^ [^Pr] = [h], neither of which is generally zero. This differs considerably 
from the constant h case treated above: in general, the equations of motion 
no longer decouple, and moreover the discrete Euler-Lagrange equations 
contain additional terms arising from the time step control function. This 
gives some insight into the delicate nature of implementing time step control 
for structure-preserving numerical integrators (as discussed at length in, e.g., 
Hairer et al. |6j). 

7. Conclusion 

7.1. Summary of results. We began this paper by developing a generalized 
theory of discrete Lagrangian mechanics, in terms of Lagrangian submani- 
folds of symplectic groupoids, which induce Poisson relations on their base 
manifolds. We also characterized the regularity and reversibility of these 
systems, providing a significant generalization of previous results. 

Applying this framework to the cotangent groupoid T*G ^ A*G of a 
groupoid G ^ Q, we were able to formulate a new theory of discrete 
constrained Lagrangian mechanics, where the discrete Lagrangian L is defined 
on a constraint submanifold N C G; this reduces to the earlier unconstrained 
theory in the special case N = G. This allows both for holonomic constraints, 
when iV is a subgroupoid of G, as well as more general nonholonomic 
constraints, and the distinction between the two is consistent with the 
continuous theory on Lie algebroids. The Lagrangian submanifold C T*G 
generated by L was shown to have the structure of an affine bundle over A^, 
associated to the conormal bundle of A^. When N is defined implicitly, by the 
vanishing of a family of functions (j)"" in a neighborhood on G, we showed that 
one can obtain natural coordinates for the fibers of S^, which correspond 
to Lagrange multipliers for the constraints. The resulting dynamics on S^,, 
therefore, specify the evolution of configurations (on the base) and Lagrange 
multipliers (on the fibers). 

In this setting, we also introduced the notion of a morphism between 
discrete constrained Lagrangian systems, which corresponds to a morphism 
of Lie groupoids that preserves the additional structure of the discrete 
Lagrangians and constraint manifolds. With this definition, it was then 
proven that morphisms of discrete constrained Lagrangian systems allow 
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for reduction of the dynamics. In addition, we studied Noether symmetries, 
and proved a discrete version of Noether's theorem, for these constrained 
systems. 

After this, we proved that, under some regularity conditions, the dynamical 
equations can also be derived from a variational principle, where the discrete 
action sum is constrained to trajectories lying in N. This establishes a 
connection between these systems and variational integrators. 

Finally, we applied the previous results to discretize several examples of 
continuous systems with constraints, including constrained mechanics and 
optimal control on Lie groups, the nonholonomic plate-ball system, and 
time-dependent systems with either fixed or adaptive time steps. 

7.2. Future directions. There are many interesting directions, regarding 
discrete constrained Lagrangian mechanics, that remain to be explored. For 
example, it is well known that if we have a regular, continuous Lagrangian 
function L : TQ — )• M, and a sufficiently small time step h > 0, then one may 
define a regular discrete Lagrangian function : Q x Q — )■ M, called the 
exact discrete Lagrangian, where (go, qi) is precisely the action integral 
along the Euler-Lagrange path from qq to qi, over a time interval of size h. 
Moreover, the exact Hamiltonian flow, for time h, is just the pushforward (via 
the discrete Legendre transformations) of the discrete flow associated with 
(cf. Marsden and West |15j). In other words, the Lagrangian submanifold 
dL^ C T* [Q X Q) is equal to the graph of the Hamiltonian flow for time h. 
This result is important for the error analysis of variational integrators on 
QxQ, and this analysis can also be extended to an exact discrete Lagrangian 
: G — )• M when the continuous Lagrangian L : AG — )• M is defined on the 
Lie algebroid of G (Marrero et al. [14j ) • It would be interesting to extend this 
construction for a constrained Lagrangian system, with constraint distribution 
A C AG, and to see if one could express such results in terms of a Lagrangian 
bisection of T*G, generated by an exact discrete Lagrangian on some N G G. 

Finally, for continuous systems, Yoshimura and Marsden [24] showed that 
a variational principle, which they call the Hamilton-Pontryagin principle, is 
closely related to Dirac structures, and thus can be quite useful for studying 
systems with nonholonomic constraints and various degeneracies. In Stern 
[15] . it was shown that a discrete Hamilton-Pontryagin principle describes the 
relationship between the generating-function and variational interpretations 
of the unconstrained discrete Lagrangian L : G — t- M. It would be interesting 
to see if this variational principle could be generalized to discrete constrained 
Lagrangians L: N C G ^ M, and how this might be connected to the 
discretization of Dirac structures and of the Courant algebroid. 
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Appendix A. Bisections of Lie groupoids 

In this appendix, we will recall some key facts about bisections of Lie 
groupoids, which provide the technical underpinning for several concepts 
appearing elsewhere in this article. Most notably, we derive an expression 
for the tangent map of Lie groupoid multiplication. This expression, which 



appeared in Xu |23j . is a crucial ingredient in the proof of the Theorem 5.1 
which we employed in the variational formulation of discrete constrained 
Lagrangian mechanics. We also apply this result to present an alternative 
derivation of the multiplication map for the cotangent groupoid. (Our 
definitions largely follow those of Cannas da Silva and Weinstein 13^ Chapter 
15].) 

A.l. Global and local bisections. Let G ^ Q he a Lie groupoid with 
source map a: G ^ Q and target map (3: G ^ Q. A submanifold 13 C G is 
called a bisection of G if the restricted maps, a|g : B ^ Q and /3|e : B ^ Q 
are both diffeomorphisms. Consequently, for any bisection B C G, there is 
a corresponding a-section Ba = (als)^^ : Q ^ G, where f3 o Ba - Q ^ Q 
is a diffeomorphism. (Mackenzie [H] defines a bisection to be the map Ba 
having this property, rather than its image B; the definitions are equivalent.) 
Likewise, there is a /3-section = (/3|b)~ : Q ^ G, where a o Bp = 
(P o Ba)~^ : Q — )• Q is a diffeomorphism. Furthermore, each bisection B C G 
defines a left action Iq: G ^ G and a right action r/g: G — )• G on the 
groupoid, given on any g £ G hy 

iisig) = Bp {a{g)) g, vsig) = gBa (/3(ff)) . 

Given two bisections 0i,02 C G, one can show that the product B1B2 = 
^Bi (^2) = 1^62 i^i) is again a bisection, and that the bisections of G in fact 
form a group. 

More generally, B C G is called a local bisection if the restricted maps 
a\B and are local diffeomorphisms onto open sets U,V C Q, respectively. 
Analogously to the global case, there exists a local a-section Ba : U ^ B 
such that (3 o Ba '■ U V is a diffeomorphism, as well as a local /3-section 
B(s:V^B such that aoBp = (j3o BaY : V ^ U is a diffeomorphism. 
(As before, an alternate but equivalent definition, cf. Mackenzie [IT], takes a 
local bisection to be the map Ba-, rather than B itself.) Each local bisection 
B C G defines a local left action £3- a~^(V) — )• a^^{U) and a local right 
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action tb'- (3^^{U) — )■ /3'^^{V), given by 

£b = Bp {a{g)) g, reig) = gB^ {P{g)) . 

It is also possible to define multiplication of local bisections; moreover, if 
i: G ^ G denotes the inversion map on G, then the inverse i{B) of a local 
bisection is again a local bisection. Thus, the local bisections of G form an 
inverse semigroup, which contains the group of global bisections as those 
elements with U = V = Q. 

A. 2. The tangent map of Lie groupoid multiplication. A key fact 
about Lie groupoids is that, for every g £ G, there exists a local bisection 
B C G such that g £ B. (In general, though, there may not exist any global 
bisection through g.) Using this property, we may now state and prove the 
following theorem on the tangent map of multiplication in G. 

Theorem A.l (Xu [23j). Let G ^ Q be a Lie groupoid, with multiplication 
mapm: G2 — )■ G. Suppose {gi,g2) G G2 is a pair of composable elements, and 
denote q = f3 (gi) = a ((72)- Then the tangent map Ti^g^ g^-^m: T(^gi.g2)G2 — ^ 
Tg^g^G is given by 

where Bi and B2 are local bisections such that gi £ Bi, g2 £ B2, and where 
= Tg^P (Vgi) = Tg^a {Vg^) £ TqQ. 

Proof. Denote by (-82)0, and (-Bi)^, respectively, the local a- and /3-sections 
corresponding to B2 and Bi. This implies that 52 = {B2)a (q) ^-^d gi = 
(9)1 so we decompose the tangent vector {vg^^Vg^) G T(^g-^ g^-jG2 as 

(9) {Vgi,Vg2) = iVg^,Tq{B2)^{Vg)) + iTg{Bl)^{Vg),Vg2) 

-{Tg{B,)^{Vg),Tg{B2)^{Vg)). 

Applying the tangent map to each of these terms individually, we observe 
for the first term that 

^{91,92)"^ (^51' (^2)^, (Vq)) = ^(9i,92)"^(^9lidG (Vgi) ,Tg, ((B2)„ o /3) (VgJ) 

= Tg, {mo{idG,iB2)^oP))ivg^). 

Now, since 

(m o (idc, (i?2)„ o 13)) (g) = m (g, (i^a), Pig)) = g (i^s)^ f3{g) = rsM, 
we conclude that 

(10) ^(91,52)"^ {Vgi^Tq {B2)^ [Vg)) = Tg^rB2 (Vgi) . 

Similarly, applying the tangent map to the second term of ([9]), it follows that 

(11) ^{9l,92)™'(^g (^l)/3 (^^g) '^^92) = ^92^Bi (^^92) • 

On the other hand, for the third term of ([o]), we have 

T(gug2MTq {Bi)p {vg) , Tg (i?2), {vg)) = T, (m o {{B^)p, (i?2)„)) {vg), 
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and since 

(m o iB2)a)) (q) = {Bi)^ (q) {B^)^ (q) = {£b, o tb, o e) (g), 

we deduce that 

(12) r(,^,,,)m(rg (Bi)^ {v,) , Tg {B2)^ K)) = Tg {iB, o rB, o e) {v,) . 



Finahy, substituting the expressions ( |10[ ), ( |11| ), (12) for the respective terms 
of Q yields the result. □ 

A. 3. Multiplication in the cotangent groupoid. Given a Lie groupoid 
G ^ Q with Lie algebroid AG — )■ Q, we have already discussed the source 
and target maps, respectively denoted a and (3, of the cotangent groupoid 
T*G ^ A*G. In this subsection, we discuss the multiplication map in, whose 
formulation depends on the properties of bisections. In particular, we derive 
an expression for this multiplication map using Theorem A.l[ 

The multiplication map fh on T*G can be characterized by the following 
two conditions (cf. Mackenzie |llj): 

(i) if ilJ'gi, lJ'g2) £ Tg^G X Tg^G is any composable pair, i.e., 13 (fig-^) = 
a ing^), then the product m (/ig^, /ig2) li^s in Tg_^g^G; 

(ii) if ivg^,Vg2) £ T^g^,g2)G2, then {m {^g^, ^ig^) ,T^g^^g^)m {vg^,Vg2)) = 

Using these properties, we can now obtain the following explicit formula for 
the cotangent multiplication map. 

Theorem A. 2. Let {fJ-gi, fJ'g2) £ ^gi^ ^ '^92^ ^ composable pair, with 
j3 {gi) = a (52) = q ^ Q o-n-d /3 (fJ-gi) = 5 (/^ga) = ^ ^qG- Define the linear 
epimorphism 

TTg : T,(^g)G ^ AgG , v^^g) ^ v,^g) - r,(g) (e o a) . 
Then the cotangent multiplication map fh is given by 

= Tg,gJi{Bi) (/^ga) + T*g^g2'^i[B2) if^9i) " (^3192 ikBi) ° ^*(i?2)) ° K) (Z^?) ' 
where Bi and B2 are local bisections of G such that gi G Bi, g2 £ B2. 
Proof. First, we calculate 

{1^92, Vg^) = {1^92, Tg2 {ii{Bl)°£Bl) (Vg^)) 

= (^9192 (^92) >^92^Bl (^^92)) 

= (^9i92^i(Si) (/^92) > ^(91,92)"^ (^91' ^92) - ^91 '•S2 (Vgi) 

+ Tg{£Bi orB^oe) (vg)), 
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where the last equahty is obtained by applying |Theorem A.l[ Observe that 

= (^92, ^91 (^i(Bi) ° rB^) (vg^) - Tg {rB^ o e) (vg)) 

= {l^92^-Te{q) {rB2 ° i) (TTg {TgJi(B^) Kl)))) 

= {^^*q^i^^g2) ,TgJi^Bi) (vg^)) , 
so altogether, we have 

(/^92'^92) = (^9*i92^i(Bl) (/^92) > ^(91,92 ^92)) 

Following essentially the same procedure, one also obtains 

(^9i'^9i) = (^gi92^i(B2) (/^9l) '^{91 ,92)"^ (^91-^^92) ) 

SO adding these together, 

(/^9l'^9i) + (M92>^^92) 

= (^9*i92^i(fii) (/^92) +^9*192 ^i(B2) (A^9i) '^(9i,92)'^(^9i'^92)> 

- (^g if^q) ,Vi(^Bi) (Vg^) + Tg^Vil^Bi) K2 ) ) • 

However, for this last term, we can use [Theorem A.l| again to write 

Vt(Bi) (Vgi) + Tg^fi^Bi) {Vgz) 

= ^9192 (^i(Bi) ° ^»(i?2)) (^9i^i?2 (^^l) + ^92^i?i (^^2)) 

= ^9192 (^i(Bi) ° ^*(i?2)) (^(91,92)"^ (f^9i'^'92)+?'g(%°'^B2°e) K)) 
= ^9192 (^i(Bi) ° ri(^B2)) (^(91,92)"^ (^^91' ^92)) + TqC {Vq) . 

But since 

{■Kg O Tge) {Vg) = Tg £ {V g) " Tg ( € O Q! O £ ) {V g) = TgE {V g) " Tg £ {V g) = 0, 

the Tgt term vanishes. Finally, we are left with 

(^9l'^9i) + (M92'^'92) = (^9*i92^i{i?i) {^^'92) + T*g^g2^i{B2) (/^9i) 

- (^9*192 ° ^i(B2)) ° A) (/^g) '^(91,92)"^ (^9i'^92)>> 



and since this equals (^m{iJ,g^, fig^) ,T(^g^g^^m{vg-^^,Vg2)), by condition (ii 
above, this completes the proof. CH 
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